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Dedicated to Claude Bardos on his 70-th birthday 
Abstract 

We consider the linear Schrodinger equation under periodic boundary 
condition, driven by a random force and damped by a quasilinear damping: 

— u + i ( — A + V(x))u = v^Au — r yn\u\ 2p u — i~/i\u\ 2q uj +\/ur](t, x). (*) 

The force rj is white in time and smooth in x. We are concerned with 
the limiting, as v — > 0, behaviour of its solutions on long time-intervals 
< t < v~ 1 T, and with behaviour of these solutions under the double 
limit t — > oo and v — > 0. We show that these two limiting behaviours 
may be described in terms of solutions for the system of effective equations 
for (*) which is a well posed semilinear stochastic heat equation with 
a non-local nonlinearity and a smooth additive noise, written in Fourier 
coefficients. The effective equations do not depend on the Hamiltonian 
part of the perturbation —i^i\u\ 2q u (but depend on the dissipative part 
— 1r\u\ 2p u). If p is an integer, they may be written explicitly. 
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Introduction 

In [KP08, KuklO] we have considered the KdV equation on a circle, perturbed 
by a random force and a viscous damping. There we suggested auxiliary effective 
equations which are well posed and describe long-time behaviour of solutions for 
the perturbed KdV through a kind of averaging. 

In this work we apply the method of [KP08, KuklO] to a weakly nonlinear 
situation when the unperturbed equation is not an integrable nonlinear PDE (e.g. 
KdV), but a linear Hamiltonian PDE with a generic spectrum. Since analytic 
properties of the latter are easier and better understood then those of the former, 
in the weakly nonlinear situation we understand better properties of the effective 
system and its relation with the original equation. Accordingly we can go further 
in analysis of long time behaviour of solutions. 

More precisely, we are concerned with z/-small dissipative stochastic pertur- 
bations of the space-periodic linear Schrodinger equation 

^-u + i(-Au + V(x))u = 0, x £ T d , (0.1) 

i.e. with equations 

^-u + iAu = uUu- lR f p {\u\ 2 )u-i ll f q (\u\ 2 )u) + vWM), xeT d , (0.2) 
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where r)(t,x) = j- t J2jLibjPj(t)ej(x). Here Au = A v u = —Au + V(x)u and the 
potential V(x) > 1 is sufficiently smooth; the real numbers p, q are non-negative, 
the functions f p (r) and f q {r) are the monomials \r\ p and \r\ q , smoothed out near 
zero, and the constants 7^, 7/ satisfy 

7i?,7/>0, 1r + 1i = 1- (0.3) 

If 1r = 0, then due to the usual difficulty with the zero-mode of the solution u, the 
term Am in the r.h.s. should be modified to A — u. The functions {ej(x), j > 1} in 
the definition of the random force form the real trigonometric base of L 2 (T d ), the 
real numbers bj decay sufficiently fast to zero when j grows, and {/3-(i), j > 1}, 
are the standard complex Wiener processes. So the noise rj is white in time and 
sufficiently smooth in x. It is convenient to pass to the slow time r = ut and 
write the equation as 

u + v~ l iAu = Au- -iRf p {\u\ 2 )u - i7// g (|w| 2 )w + i](t, x), (0.4) 

where ii = du/dr. The equation is supplemented with the initial condition 

u(0,x) = u {x). (0.5) 

It is known that under certain restrictions on p, q and d the problem (0.4), (0.5) 
has a unique solution u u (t,x), t > 0, and eq. (0.4) has a unique stationary 
measure [i v . We review these results in Section 1 (there attention is given to the 
Id case, while higher- dimensional equations are only briefly discussed). 

Let {ipk, k > 1}, and k > 1}, be the eigenfunctions and eigenvalues of 

Ay, 1 < Ai < A2 < We say that a potential V is nonresonant if J2^jLi ^j s j 

for every finite non-zero integer vector (si, ■ ■ ■ )• In Sections 1.4, 1.5 we show 
that nonresonant potentials are typical both in the sense of Baire and in the sense 
of measure. Assuming that V is nonresonant we are interested in two questions: 

Ql. What is the limiting behaviour as v — > of solutions u y (r,x) on long time- 
intervals < r < T? 

Q2. What is the limiting behaviour of the stationary measure [i v as v — > 0? 

For any complex function u(x), x G T d , denote by ^/(u) = v = (i>i,i>2, • • • ) 
the complex vector of its Fourier coefficients with respect to the basis {(fk}, i.e. 
u(x) = J^Vjipj. Denote 

Ij = ^\vj\ 2 , <fj = AigVj, j > 1. (0.6) 
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Then (I, ip) E x T°° are the action-angles for the linear equation (0.1). The 
v- and (/, <y?)-variables are convenient to study the two questions above. Writing 
(0.4) in the (J, y?)-variables we arrive at the following system: 

£,, = ..., = + (0.7) 

where the dots stand for terms of order one (stochastic and deterministic). We 
have got slow/fast stochastic equations to which the principle of averaging is for- 
mally applicable (e.g., see [AKN89, LM88] for the classical deterministic averaging 
and [Kha68, FW98] for the stochastic averaging). Denoting ij(r) = Ij{u u {r)) and 
averaging in the /-equations in (0.7), using the rules of the stochastic calculus 
[Kha68, FW98] and following the arguments in [KP08], we show in Section 2 that 
along sequences Uj — > we have the convergences 

V(I^(.))^V(I%)), (0.8) 

where the limiting process I°(t), < r < T, is a weak solution of the averaged 
/-equations. As in the KdV-case the averaged equations are singular and we do 
not know if their solution is unique. So we do not know if the convergence (0.8) 
holds as v — > 0. To continue the analysis we write eq. (0.4) in the v- variables 

v k + iv~ x \ k v k = P k (v) + B kjPj(r), (0.9) 

where the drift P k and the dispersion Bkj are written explicitly in terms of the 
r.h.s. of eq. (0.4). It turns out that the Hamiltonian term — i f q {\u\ 2 )u con- 
tributes to P{v) a term which disappears in the averaged /-equations. We re- 
move it from P(v) and denote the rest P(v). For any vector 9 = (9i,9 2 ,...) G T°° 
denote by the linear transformation of the space of complex vectors v which 
multiplies each component Vj by e %di . Following [KuklO] we average the vector 
field P by actions of the transformations and get the effective drift R(v) = 
f JOO $_6)P($6if ) d9. In Section 3.1 we show that 

R k (v) = -\ k v k + R° k (v), (0.10) 

where R°(v) is a smooth locally Lipschitz nonlinearity. 

Since the noise in (0.9) is additive (i.e., the matrix B is v- independent), then 
construction of the effective dispersion, given in [KuklO] for non-additive noise, 
simplifies significantly and defines the effective noise for eq. (0.9) whose k-th 
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component equals tf^li) d/3 k (r). Accordingly the effective equations for 
(0.4) become 

v k = R k (v)dr+ (j2b^ii) 1/2 d(3 k (T), k > 1. (0.11) 
i 

By construction this system is invariant under rotations: if v(t) is its weak solu- 
tion, then <&qv{t) also is a weak solution. Due to (0.10) this is the heat equation 
u = —Au for a complex function u(t,x), perturbed by a non-local smooth non- 
linearity and a non-degenerate smooth noise, written in terms of the complex 
Fourier coefficients Vj. It turns out to be a monotone equation, so its solution is 
unique (see Section 3.2). 

In particular, if in (0.4) p = 1, then the system of effective equations takes 
the form 

/ \ 1/2 

v k = -v k ((X k -M k )+ lR J2N 2 Lki)dr+( K J2 b i^li) d Pk(r), k > 1, (0.12) 

i 

where M k = f V(x)(p 2 (x) dx and L kl = (2 — 5 k i) f ipl(x)<p 2 (x) dx. See Exam- 
ple 3.1 (the calculations, made there for d — 1, remain the same for d > 2). 

It follows directly from the construction of effective equations that actions 
{I(v k (r)) = ^\v k (r)\ 2 , k > 1} of any solution v(t) of (0.11) is a solution of the 
system of averaged /-equations. On the contrary, every solution I°(r) of the 
averaged /-equations, obtained as a limit (0.8), can be lifted to a weak solution 
of (0.11). Using the uniqueness we get 

Theorem 0.1. Let I v {r) = I(u v (t)), where u v (r), < r < T ', is a solution of 
(0.4), (0.5). Then\im u ^ V(I u (-)) = V(I°(-)), where I°(r), < r < T, is a weak 
solution of the averaged I-equations. Moreover, there exists a unique solution 
v(t) of (0.11) such that v(0) = v = #(u ) and V(I(v{-)) = V(I°(-)), where 

lHr))^\\vAr)\ 2 - 

The solutions /°(r) and v (r) satisfy some apriori estimates, see Theorem 3.5. 
Concerning distribution of the angles (p{u v (r)) and their joint distribution with 
the actions see Section 2.4. 

Now let \i v be the unique stationary measure for eq. (0.4) and u' v be a corre- 
sponding stationary solution, V{u' u {r)) = /i u . As above, along sequences Uj — > 
the actions I'" 3 (r) = I{u' V3 {t)) converge in distribution to stationary solutions 
I'{t) of the averaged /-equations. These solutions can be lifted to stationary 
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weak solutions v'{t) of effective equations (O.lf ). Since that system is monotone, 
then its stationary measure m is unique. So the limit above holds as v — > 0. As 
the effective system is rotation invariant, then in the (/, </?)- variables its unique 
stationary measure has the form dm = mi(dl) x dip. It turns out that the mea- 
sure lim„^ fJ>" also has the rotation-invariant form and we arrive at the following 
result (see Theorem 4.3, 4.4 for a precise statement): 

Theorem 0.2. When u — >■ we have the convergences V (I{u' h '(•)) — 1 VI{v'{-)) 
and \& o n v -± m, where dm = mj(dl) x dip. 

Accordingly every solution u"(t) of (0.2) obeys the following double limit 

lim lim V(u v (t)) = o m. (0.13) 

By Theorems 0.1 and 0.2, the actions I(u u (t)) of a solution u u of (0.4), (0.5) 
converge in distribution to those of a solution v(t) of the effective system (0.11) 
with v (0) = ty(uo), both for < r < T and when r — > oo. We conjecture that 
this convergence hold for each r > 0, uniformly in r (the space of measures is 
equipped with the Wasserstein distance). 

In Example 4.6 we discuss Theorem 0.2 for equations with p — 1, when 
the effective equations become (0.12). In particular, we show that Theorem 0.2 
implies that in equations (0.2) with small v there is no direct or inverse cascade 
of energy. 

In Example 4.5 we discuss Theorem 0.2 for the case 7_r = (when the non- 
linear part of the perturbation is Hamiltonian) and its relation to the theory of 
weak turbulence. 

We note that the effective equations (0.11) depend on the potential V(x) in 
a regular way and are well defined without assuming that V(x) is nonresonant 
(cf. equations (0.12)). In particular, if V M (x) — > 1 as M — > oo, where each 
V M (x) > 1 is a non-resonant potential, then in (0.13) m M — m(l), where m(l) 
is a unique stationary measure for eq. (0.12) with V{x) = 1. In this equation 
^ki = 8k,h M k = 1 and the constants L k[ can be written down explicitly. 

In Section 5 we show that Theorems 0.1, 0.2 remain true for Id equations with 
non- viscous damping (when Aw in the l.h.s. of (0.2) is removed, but 7r > 0). 

Inviscid Limit. A stationary measure \x v for eq. (0.4) also is stationary for 
the fast-time equation (0.2). Let U"(t) be a corresponding stationary solution, 
VU u (t) = fi u . It is not hard to see that the system of solutions U u (t) is tight on 
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any finite time-interval [0,T]. Let \JJ v \Vj — > 0}, be a converging subsequence, 
i.e. 

V(U V ') -± Q*, fi Vj -± n*. 

Then fj,* is an invariant measure for the linear equation (0.1) and Q* = D(U *(•)), 
where U*(t),0 < t < T, is a stationary process such that / D(U*(t)) = //* and 
every trajectory of U* is a solution of (0.1). The limit V(U l/j ) V(U*) is 
the inviscid limit for eq. (0.2). Eq. (0.1) has plenty of invariant measures: if 
we write it in the action-angle variables (0.6), then every measure of the form 
m(dl) x dip is invariant (see [KS04] for the more complicated inviscid limit for 
nonlinear Schrodinger equation). Theorem 0.2 explains which one is chosen by 
eq. (0.2) for the limit lim^ (/,". 

The inviscid limit for the damped/driven KdV equation, studied in [KP08, 
KuklO] is similar: the limit of the stationary measures for the perturbed equa- 
tions is a stationary measure of the corresponding effective equations. Due to a 
complicated structure of the nonlinear Fourier transform which integrates KdV, 
uniqueness of their invariant measure is not proved yet. So the final results con- 
cerning the damped/driven KdV are less complete than those for the weakly 
perturbed CGL equation in Theorem 0.2. 

Finally consider the damped/driven 2d Navier-Stokes equations with a small 
viscosity v and a random force, similar to the forces above and proportional to 

v' t -uAv + {v-V)v + Vp = v ^ri{t,x); divv = 0, v G R 2 , xeT 2 (0.14) 

It is known that (0.14) has a unique stationary measure //", the family of measures 
{//", < v < 1} is tight, and every limiting measure lim^._^. /i^ is a non-trivial 
invariant measure for the 2d Euler equation (0.14)^ =0 , see Section 5.2 of [KS10]. 
Hovewer it is non-clear if the limiting measure is unique and how to single it out 
among all invariant measures of the Euler equation. The research [KP08, KuklO] 
was motivated by the belief that the damped/driven KdV is a model for (0.14). 
Unfortunately, we still do not know up to what extend the description of the 
inviscid limit for the damped/driven KdV and for weakly nonlinear CGL in terms 
of the effective equations is relevant for the inviscid limit of the 2d hydrodynamics. 

Agreements. Analyticity of maps B\ — > B2 between Banach spaces B\ and B2, 
which are the real parts of complex spaces B\ and B^, is understood in the 
sense of Frechet. All analytic maps which we consider possess the following 
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additional property: for any R a map analytically extends to a complex (5r > 0)- 
neighbourhood of the ball {H^ < R} in B\. 

Notations, xa stands for the indicator function of a set A (equal 1 in A and equal 
outside A). By x(t) we denote various functions of t such that x(t) — > when 
t — >■ oo, and by Xoo(t) denote functions x(i) such that x(i) = o(t~ N ) for each N. 
We write x(t) = x(t; R) to indicate that x(t) depends on a parameter R. 

Acknowledgments. I wish to thank for discussions and advice Patrick Gerard, 
Sergey Nazarenko, Andrey Piatnitski and Vladimir Zeitlin. This work was sup- 
ported by l'Agence Nationale de la Recherche through the grant ANR-10-BLAN 
0102. 



1 Preliminaries 

1.1 Apriori estimates. 

We consider the Id CGL equation on a segment [0, 7r] with a conservative linear 
part of order one and a small nonlinearity. The equation is supplemented with 
Dirichlet boundary conditions which we interpret as odd 27r-periodic boundary 
conditions. Introducing the slow time r = vt (cf. Introduction) we write the 
equation as follows: 

u + iv~ x ( — u xx + V(x)uj = xu xx — 7 R |w| 2r u — i^j\u\ 2q u 

d (1 1) 

+ — ^ bj(3j(r)ej(x), u{x) = u(x + 2-k) = —u{—x). 

i=l 

Here u = -^u, p,qEZ + := NU {0} (only for simplicity, see next section), x > 0, 
constants 7^ and 7/ satisfy (0.3) and R 3 V(x) > is a sufficiently smooth even 
27r-periodic function, {ej,j > 1} is the sine-basis, 

eAx) = —= sin jx, 
v 7r 

and f3j,j > 1, are standard independent complex Wiener processes. That is, 
(3j(r) = f3j(r) + i/3-j(r), where P±j{r) are standard independent real Wiener 
processes. Finally, the real numbers bj all are non-zero and decay when j grows 
in such a way that Bi < 00, where 

00 

B r := 2 f rh ) < 00 for r > °- 
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By W, reKwe denote the Sobolev space of order r of complex odd periodic 
functions and provide it with the homogeneous norm || • || r , 



\u\\l 



EN 2p ' Hlo 



1=1 



(ifreN, then N| r = |f£| La ). 

Let w(t, x) be a solution of (1.1) such that u(0, x) = uq. Applying Ito's formula 
to |||m|| 2 we get that 

d QlMI 2 ) = (-Tr-kHS - *||? + ^ )dr + dM(r), (1.2) 

where M(r) is the martingale f^YlbjUj ■ d/3j(r). Here |u| r stands for the L r - 
norm, 1 < r < oo, and for complex numbers z±, z 2 we denote by z\ • z 2 their real 
scalar product, 

z± • z 2 — Re z{zi. 

So (uj + iu-j) ■ (dj3j + id/3-j) = Ujd(3j + u_jdf3_j. From (1.2) we get in the usual 
way (e.g., see Section 2.2.3 in [KS10]) that 

Ee p»\HT)r < c ^ 5o) || Mo ||) Vr > o (i.3) 

for a suitable p x > 0, uniformly in v > 0. 
Denoting 

£{r) = \\\u{r) || 2 + 7r jT \u\XXlds + |jT \\u\\lds 

and noting that the characteristic of the martingale M is (M)(r) = J2tf\uj\ 2 < 
b 2 M \\u\\ 2 , where bu = max \bj\, we get from (1.2) that 

£{t) < \\\u \\ 2 + \B t + M(t) - f \\u\\ 2 ds 

2 J 



< lll«ol| 2 + hB r + H^b 



-U2 
M 



(x6 M 2 M(r))-i(x6 M 2 M)(r) 



Applying in a standard way the exponential supermartingale estimate to the term 
in the square bracket in the r.h.s. (e.g., see [KS10], Section 2.2.3 ), we get that 

P{sup(£(r) - \B q t) > ±|M| 2 + p} < e~ 2 ^, (1.4) 

r>0 
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for any p > 0. 

Now let us re- write eq. (1.1) as follows: 

u + ii/~ 1 (-u xx + V(x)u + i>7 I \u\ 2q u) = xu xx - 7r\u\ 2p u+-^- bjPj(T)ej. (1.5) 
The l.h.s. is a Hamiltonian system with the hamiltonian —h>~ l H{u), 

H(u)= 1 -(Au,u)+ ll ^ J \u\ 2q + 2 dx, A=-^ + V(x). 
For any j G N we denote 

\\u\\' r 2 = (A r u, u). 
Then dH(u)(v) = (Au,v) + 7/z/(|w| 2 %, v) and 

oo 

5'2^ b 2 d 2 H{u){e v e,) = IfiJ + 7/ z/X(r), 

where 



2 



and 

X{r) = 2pRe j (\u\ 2q - 2 u 2 J2b 2 e j (x) 2 )dx + j \u\ 2q ^b 2 ej{x) 2 dx 



<CB \u(t)Q. 



Therefore applying Ito's formula we get that 

dH{u{r))= {{- lR (AuM 2p u) + x(Au,u xx ) - 7j i/ 7fl | | M | 2 ^ +2 rfx 
+ x7/z/(|m| 2<? m, m^) + Si + 7/z/X(r)jrfr + dM(r), 

where dM(r) = Yl bj(Au + 7/z/|w| 2 %, e,) • d(3Ar). 

Denoting U q (x) = -^r[U q+1 and U p (x) = -^u p+l , we have 

(\u\ 2q u,u xx ) <- J \u x \ 2 \u\ 2q dx = -ll^^ll 2 , 



(1.6) 
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and a similar relation holds for q replaced by p. Accordingly, 

dH{u{r)) < --{x\\u\\l + 7r||^ p || 2 + x 7/ z/||— f/ g || 2 

+viil R J \u\ 2p+2q+2 dx - C„\\u\\ 2 - 2BCj dr + dM(r), 

where C x may be chosen independent from x if 7 _r > 0. Considering relations 
on H{u) m , m > 1, which follow from (1.7) and (1.6), using (1.4) and arguing by 
induction we get that 

e( sup H(u(t)) m + ^ [ T H^WuWlds] < H(u ) m +C m (x,T,B 1 )(l+\\u \\ c ™), 

\0<t<T 2. J J 

(1.8) 

EH(u(t)) m < a ra (^,5 1 )(l + iJ(M ) m + ||M || Cm ) Vt>0. (1.9) 

Estimates (1.8) in a traditional way (cf. [Hai02, KS04, Oda06, Shi06]) imply that 
eq. (1.1) is regular in space "H 1 in the sense that for any u G "H 1 it has a unique 
strong solution, satisfying (1.4), (1.8) 



1.2 Stationary measures. 

The a-priori estimates on solutions of (1.1) and the Bogolyubov-Krylov argument 
(e.g., see in [KS10]) imply that eq. (1.1) has a stationary measure fi u , supported 
by space H 2 . Now assume that 

bj + Vj. (1.10) 

Then the approaches, developed in the last decade to study the 2d stochastic 
Navier-Stokes equations, apply to (1.1) and allow to prove that under certain 
restrictions on the equation the stationary measure fi u is unique. In particular 
this is true if 7/ = (the easiest case), or if p > q and 7_r 7^ (see [Oda06]), or 
if 1r — and p = 1 (see [Shi06]). In this case any solution u{t) of (1.1) with 
u(0) = uo G "H 1 satisfies 

Vu{t)^ii v as t^oo. (1.11) 
This convergence and (1.3), (1.9) imply that 

J e pM ^{du) < C{h,B), (1.12) 
J |M| 2 "V(dw) < C m (x,Si) Vm. (1.13) 
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1.3 Multidimensional case. 

In this section we briefly discuss a multidimensional analogy of eq. (1.1): 
ii + iv~ l Au = Au - i R fp{\u\ 2 )u - ?7//g(|w| 2 )w 

+ ^f2 b j(3 j (r)e j {x), u = u{t,x), xeT d . (L14) 

Here Au = -Au + V(x)u, V G C JV (T i ,R) and V(x) > 1. The numbers 77, 7 R 
satisfy (0.3). Functions f p > and f q > are real-valued smooth and 

/ p (0 = t p for t > 1, = for * > 1, 

where p, g > 0. If 7^ = 0, then the term Au in the r.h.s. should be modified 
to A — u. By {e k , k > 1}, we denote the usual trigonometric basis of the space 
L 2 (T d ) (formed by all functions ir^ d ^ 2 f Sl (07) . . . f Sd (x d ), where each f s (x) is sin sx 
or cossrr), parameterised by natural numbers. These are eigen-functions of the 
Laplacian, — Ae r = \ r e r . We assume that 

B' Ni =2j2^bl<oo, (1.15) 

k 

where Ni = Ni(d) is sufficiently large. In this section we denote by (W, || ■ || r ) 
the Sobolev space W = H r (T d ,C), regarded as a real Hilbert space, and (•, •) 
stands for the real L 2 -scalar product. 

Noting that (f p (\u\ 2 )u — \u\ 2p u) and (f q (\u\ 2 )u — \u\ 2q u) are bounded Lipschitz 
functions with compact support we immediately see that the a-priori estimates 
from Section 1.1 remain true for solutions of (1.14). Accordingly, for any u G 
"H 1 nL 2q+ 2 eq. (1.1) has a solution u(t,x) such that u(0, x) = u , satisfying (1.3), 
(1.8), (1.9). 

Now assume that 

2 

p,q< 00 if d= 1,2, p,q<- if d > 3. (1.16) 

d 2 

Applying Ito's formula to the processes (A m u(T), u(r)) n , m,n > 1, using (1.3), 
(1.8), (1.9) and arguing by induction (first in n and next in m) we get that 

E( sup ||«(r)||^+ / T ||^)irL +1 ||^)|lir^, n 17 , 

\0<t<T Jo J (1-17) 

< \Wo\\'ll + C(m,n,T)(l + \\u \\ c ™'") : 
12 



E||«(r)|| 



/ 2n 
2m 



<C(m,n) Vr>0, 



(1.18) 



for each m and n, where C(m,n,T) and C(m,n) also depends on |V| c jv and B Nl 
(see (1.15)), and N = N(m), Ni = Ni(m). 

Relations (1.17) with m = mo > 1 in the usual way (cf. [Hai02, KS04, Oda06, 
Shi06]) imply that eq. (1.14) is regular in the space T-L m ° C\L2 q +2 in the sense that 
for any u G H m ° ni 29 +2 it has a unique strong solution u(t, x), equal u at t = 0, 
and satisfying estimates (1.3), (1.17) with m = m for any n. By the Bogolyubov- 
Krylov argument this equation has a stationary measure supported by the 
space % m ° PI L 2?+ 2, and a corresponding stationary solution u u {r), Vu v {r) = ji v , 
also satisfies (1.3) and (1.18) with m = mo. 

If (1.10) holds and (1.16) is replaced by a stronger assumption, then a sta- 
tionary measure is unique. If 7/ = 0, the uniqueness readily follows, for example, 
from the abstract theorem in [KS10]. In [Shi06] this assertion is proved if 

7^ = and q< 1 if d= 1, q < 1 if d = 2, q < 2/d if d = 3. (1.19) 

In [Oda06] it is established if 



the argument of this work also applies if p > q. 

Note that when 7^ = or when p < q (i.e., when the nonlinear damping is 
weaker than the conservative term), the assumptions (1.19), (1.20), needed for the 
uniqueness of the stationary measure, are much stronger than the assumptions 
(1.16), needed for the regularity. This gap does not exist (at least it shrinks a 
lot) if the random force in eq. (1.14) is not white in time, but is a kick-force. See 
in [KS00] the abstract theorem and its application to the CGL equations. 

1.4 Spectral properties of Ay\ one-dimensional case. 

As in Section 1.1 we denote Ay = A = —d 2 /dx 2 + V(x), where the potential 
V(x) > belongs to the space of C^-smooth even and 27r-periodic functions, 
N > 1. Let ipi, (f2, ■ ■ ■ be the L 2 -normalised complete system of eigenfunctions 
of A v with the eigenvalues 1 < Ai < A 2 < Consider the linear mapping 



P = Q, In, li > if d = 1, 2, and p = q < 



2 



, 7fl , 7/ >0 if d>3; (1.20) 



d-2 



* : U 3 u{x) h-> v = (v u v 2 , . . . ) e C°°, 
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defined by the relation u(x) = Yl v k i Pk{ x )- in the space of complex sequences v 
we introduce the norms 

k>i 

and denote h m = {v \ \v\hm < 00} . Due to the Parseval identity, \& : H — > h° is 
a unitary isomorphism. For any m G N we have |t>|| m = (A m u(x), u(x)). So the 
norms |^|^ m and II XL 1 1 fn equivalent for m — 0, . . . , N. Since ^* = * 1 , then the 
norms are equivalent for integer \m\ < N. By interpolation they are equivalent 
for all real \m\ < N. So 

the maps ^ : V, m — >■ /i m , |m| < TV, are isomorphisms. (1-21) 

Denote G = : /i m ->■ H m . Then 

$oioG = diag{A fe , k > 1} =: A 

Consider the operator 

£ := o (-A) oG = f o(i-l/)oG = l-f oyoG=:A-£°. (1.22) 

By (1.21) £° = ^ o V o G defines bounded maps 

C°:h m ^h m y\m\<N, (1.23) 

and in the space h° it is selfadjoint. 

For any finite M consider the mapping 

A M :C?^R M , V(x)^(X 1 ,...,X M ). 

Since the eigenvalues Xj are different, the mapping is analytic. As the functions 
<p\, ipl, . . . are linearly independent by the classical result of G. Borg, then for 
any V G the linear mapping 

dA M (V) : Cf -> R M is surjective. (1.24) 

In the space consider a Gaussian measure hk with a non-degenerate corre- 
lation operator K (so for the quadratic function f(V) = (V,£,)l 2 (V,v)l 2 we have 
/ f(V)fx K (dV) = (K£,ri)). Relation (1.24) easily implies 



14 



Lemma 1.1. For any M > 1 the measure A M o ji K is absolutely continuous with 
respect to the Lebesgue measure on M M . 

We will call a vector A e M°° nonresonant if for any non-zero integer vector s 
of finite length we have 

A-s^O. (1.25) 

A potential V(x) is called nonresonant if its spectrum A(V) = (Ai,A 2 ,...) is 
nonresonant. The nonresonant potentials are defined in by a countable family 
of open dense relations (1.25). So 



the nonresonant potentials form a subset of C, 
of the second Baire category. 



N 

e 



;i.26) 



Applying Lemma 1.1 we also get 

the nonresonant potentials form a subset of of full /Ik measure, (1-27) 

for any Gaussian measure /Ik as above. 

The non-resonant vectors A are important because of the following version of 
the Kronecker-Weyl theorem: 

Lemma 1.2. Let f e C n+1 (Y n ) for some n e N. Then for any nonresonant 
vector A we have 

hm 1 [ T f(q + tA n )dt=(2n)- n [ f dx, A n = (A 1 , . . . , A„), 

uniformly in qo aY n . The rate of convergence depends on n, A and |/|c n+1 - 

Proof. Let us write f(q) as the Fourier series f(q) = Ylfs els q - Then for each 
non-zero s we have |/ s | < C n+ i|/|cn+i|s| _n_1 . So for any e > we may find 



R = R F such that 



\s\>RJs t 



<\s\>R 
cT 



< | for each q. Now it suffices to show that 



f R (q + tA n )dt-f 



~ 2 



VT > T F 



;i.28) 



for a suitable T £ , where / fi (g) = J2\ 8 \<r fs£ %s ' q - But 
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< 



TIs- A n |' 



for each nonzero s. Therefor the l.h.s. of (1-28) is 

< \ Qj^l*- A»|) '^l/.l <T^\j\ c ,C{R,tC). 
Now the assertion follows. □ 

1.5 Spectral properties of Ay\ multi-dimensional case. 

Now let, as in Section 1.3, A = A v be the operator A = —A+V(x), x e T d , where 
1 < V(x) E C N (T d ). Let {<f k (x), k > 1} be its L 2 -normalised eigenf unctions and 

{Afc, k > 1}, be the corresponding eigenvalues, 1 < Ai < A 2 < For any 

M > 1 denote by F M C C N (T d ) the open domain 

i^M = {V | Ax < A 2 < ••• < A M }. 

Its complement F M is a real analytic variety in C N (T d ) of codimension > 2, so 
F M is connected (see [KK95] and references therein). The functions Ai, . . . , \m 
are analytic in F M . Let us fix any non-zero vector s G Z°° such that S/ = for 

1 > M. The set 

Q s = {V e F M | A(F) -s = 0} 

clearly is closed in Fm- Since the function A(V) ■ s is analytic in Fm, then either 
Q s = Fm, or Q s is nowhere dense in Fm- Theorem 1 from [KK95] immediately 
implies that Q s ^ Fm, so (1.26) also holds true in the case we consider now. 

Let hk be a Gaussian measure with a non-degenerate correlation operator, 
supported by the space C N (T d ). As A(V) ■ s is a non-trivial analytic function on 
Fm and F M is an analytic variety of positive codimension, then (j>k(Qs) = (e.g., 
see Theorem 1.6 in [AKSS07]). Since this is true for any M and any s as above, 
then the assertion (1.27) also is true. 

2 Averaging theorem. 

The approach and the results of this section apply both to equations (1.1) and 
(1.14). We present it for eq. (1.1) and at Subsection 2.5 discuss small changes, 
needed to treat (1.14). Everywhere below T is an arbitrary fixed positive number. 
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2.1 Preliminaries. 

In eq. (1.1) with u G V, 1 we pass to the f-variables, v = ty(u) G h 1 : 

v k + iv- 1 \ k v k = Pk^dT + ^Bkjdfijir), k > 1. ^.i) 

Here -B/y = ^fcj&j (a matrix with real entries, operating on complex vectors), and 

P k = Pi + P 2 k + Pi (2.2) 

where P 1 , P 2 and P 3 are, correspondingly, the linear, dissipative and Hamiltonian 
parts of the perturbation: 

where w = G(t>). We will refer to equations (2.1) as to the v-equations. 

For k > 1 let us denote = /(%) = ||^fc| 2 and <p k = <p(v k ) =Argv k G S 1 , 
where ip(0) = G S 1 . Consider the mappings 

Uj:h r 3v^I = (h,I 2 ,...)e h r I+ , Tl^:h r 3v^ip = <p 2 , . . . ) G T°°. 
Here h r I+ is the positive octant in the space 

h r i = {I\\I\ h r=2j2j 2 V 1 \<00}. 

3 

We will write 

nj(*(«)) = /(«), n v (*( u )) = ^(u), (n 7 x n^)(*(«)) = (/ x v )(u). 

The mapping / : W — > h r j is 2-homogeneous continuous, while the mappings 
(f :W — >■ T°° and (Ixip) :W^-h r I x T°° are Borel-measurable and discontinuous 
(the torus T°° is given the Tikhonov topology and a corresponding distance). 

Now let us pass in eq. (2.1) from the complex variables v k to the real variables 
h > 0, <p k G S 1 : 



d/fc(r) = K-P fc )(^)rfr + F fc 2 rfr + ^^(^-ciA), >fc = Vz] 6 ^' (2.3) 
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and 

dVk(r) = (V^fc + \v k \~ 2 (iv k ) ■ P k - \v k \~ 2 ^ b i(^ki • Ufc)(*w ■ iv k )j dr 

i 

+ ^2\vk\~ 2 k^ k i(iv k - d(3 t ) / 24 x 
i 

=: (i/^Afc + G fc (t;)) rfr + J>^) (^j • d/3,(r)) . 

Due to (1.22), (1.23) 

P(y) = xAv + P°(y), P° : h r ^ h r V § < r < iV, 

where the map P° is real analytic. The mapping P°(v) and its differential dP°(v) 
both have a polynomial growth in \v\hr. Therefore 

\P(v) - P(u m )lr-2-i/3 < m-^Qdul^), 

where Q is a polynomial. Here for any v = (i>i,i>2, ■ ■ ■ ) G /i 1 we denote v m = 
(vi, . . . , v m ) G C m and identify it with the vector (v\, . . . .v m , 0, . . . ) G h 1 . 

The functions G k and g k i are singular as v k = and satisfy the following 
estimates: 

\ G k(v)X{\v k \>6}\ < S^QMr), (2.5) 
\9kj(v)X{\v k \>8}\ < S~ l r N QkN{\v | r ), (2.6) 

where Q/e and Qfciv are polynomials. 

For any vector 6 = (6i,6 2 , . . .) £ T°° we denote by the unitary rotation 

$ e : /i r ->■ h r , v Vg, where w 0j = e t6j Vj V j. 

By (F) etc we denote the averaged functions, (F)(v) = f Jao F(<&gv) d6. They are 
(^-independent, so (F) = (F)(Hi(v)). The functions (P), (F), . . . also satisfy the 
estimates above. So 

\((v k ■ P k ))(I m ) - ((v k ■ P k ))(I)\ < m^ 3 C k Q(\I\ h] ), 

where Q is a polynomial. 

Since the dispersion matrix {B k j} is non-degenerate, then repeating for equa- 
tions (2.1) and (2.3) the arguments from Section 7 in [KP08] (also see Section 6.2 
in [KuklO]), we get 
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Lemma 2.1. Let v u (t) be a solution of (2.1) and I v {r) = I(v"(t)). Then for 
any k > 1 the following convergence hold uniformly in v > 0: 

[ P{I k (r) <5}dr^0 as 5 ->■ 0. (2.7) 
Jo 

(Certainly the rate of the convergence depends on k.) 

2.2 The theorem. 

Let us abbreviate 

h} = h, h 1 I = h I , C([0,T\,h I+ )=Hi, 

where hi + is the positive octant {I E hj \ Ij > Vj}. Fix any u e h. Due to 
the estimates (1.1) and (1.9) the set of laws {V{I U {■))}, < v < 1, is tight in 
Denote by Q° any limiting measure as z/ = i/j — > 0, i.e. 

p(/^(0)-"Q° as ^-->o. 

Let us consider the averaged drift (((ffc • Pk))(I) + ^ 2 ) for eq. (2.3). We 
have 

((v k -Pk))(v)= f (e idk v k )-P k (<£> e v)d9 = v k -R k (v), (2.8) 

where i?^, = f TOO (e~ l0k P k (Qgv)) d6 (note that ((v k • P k )) depends only on I = 
11/(1;), while R' k (v) depends on v). The diffusion matrix for (2.3) is {A kr , k,r > 
1}, where 

Its average is 

(A fcr )( V ) = J2 b * [ Re (e^- 9 ^)^ 
^ </t°° 



1/2 (2.9) 



|2 



= <Jfcr|Ufc| 2 lfc 2 , = \ J2tf\y kl \ 

Due to (1.21), 

J2 Y kk 2m <C m B m Vm. (2.10) 
Our first goal is to prove the following averaging theorem: 
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Theorem 2.2. The measure Q° is a solution of the martingale problem in the 
space hj with the drift ((vk-Pk){F)+Y k )dr and the diffusion matrix (A kr )(I). That 
is, Q° = T>(I°(-)), where the process I°(t) is a weak (in the sense of stochastic 
analysis) solution of the averaged equations 

dh = «K • P k )}(I) + Y 2 ) dr + J2 (Vifyjl) d ^(T), k > 1; (2 . n) 

r 

1(0) = I = Ui(v ). Moreover, 

E sup |/ (r)IL <Cn(KHf + l) Vn, (2.12) 



0<t<T 



E 



[ T \I (r)\ hj dr<C(\\u \\l + l). (2.13) 
Jo 

Proof. The crucial step of the proof is to establish the following lemma: 

Lemma 2.3. Let F(v) be an analytic function on the space h = h 1 which extends 
to an analytic function on h 2 ^ 3 of a polynomial growth. Then 



2T := E max 

0<r<T 



f (F(r( s ),^(s))-(F)(r(s)))ds 

Jo 



->• as v ->■ 0. (2.14) 



The lemma is proved below in Section 2.3, following the argument in [KP08]. 
Now we derive from it the theorem. Let us equip the space Hi with the Borel 
sigma-algebra F, the natural filtration of sigma- algebras {F T1 < r < T} and 
the probability Q°. The fact that the processes I%(t) - £ ({v v k ■ P k ){v v ) + Y k 2 ) ds 
are martingales, the convergence V(I Uj (-)) — ^ Q° and Lemma 2.3 with F = F k 
imply that the processes Z k {r) = Ikij) — Jq ((( v k • Pk))(H s ) + ^k) k > 1, 
are <5°-martingales, cf. Section 6 of [KP08]. 

Similar to (2.14) we find that 



E max 

0<t<T 



(F(F(s),^(s))-(F)(F(s)))ds 



->■ as v ->■ 0. 



Then using the same arguments as before, we see that the processes Zk(r)Zj(r) — 
J {Akj){I{s)) ds also are Q°-martingales. That is, Q° is a solution of the martin- 
gale problem with the drift (F k ) + Y k and the diffusion (A). Hence, Q° is a law 
of a weak solution of eq. (2.11). See [Yor74]. 

Estimates (2.12), (2.13) follow from (1.8) and the basic properties of the weak 
convergence since ~ \v |^ = Hi(v)\h m - □ 
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2.3 Proof of Lemma 2.3. 

Fix any m > 1 and denote by I u ' m , (p u > m etc the vectors, formed by the first m 
components of the infinite vectors F, tp v , etc. Below R denotes a suitable function 
of v such that R(u) — > oo as v — > 0, but 

i^z/^O as z/-»0, VAT. (2.15) 

Denote by Qr = Q R the event 

n R = { sup y(r)\ hl >R}. 

0<t<T 

Then P(f2^) < x 00 (R) uniformly in v (see Notations). We denote 

Pn R (Q) = P(^ C R n Q), E QR (f) = E(f X a k ). 

Since for \v\ h i < R we have |t> — v m \ h 2/3 < C(R)m~ 1 / ? ' and since F is Lipschitz 
on /i 2 / 3 , then 



21" < x oc (i?) + C k (R)m- 2/3 + E Qr max 

0<r<T 



/ (P(r> m , ip v ' m ) - (F) m (r' m )) ds 

Jo 



Here (F) m stands for averaging of the function T m 3 I m h-» F(I m , 0, . . . ). So it 
remains to estimate for any m and i? an analogy 21^ R of the quantity W for the 
finite-dimensional process F> m (j) on the event £l R (where its norm is < R), 

Wm,R = E« fi max f (P(F ,m , - (F) m (J^)) ds . 

Consider a partition of [0, T] by the points 

Tj=To+jL, 0<j<K, 
where tk is the last point Tj in [0, T). The diameter L of the partition is 

L = \fv, 

and the non-random phase r G (0, L] will be chosen later. Denoting 

rjt = f (F(r> m , tp u > m ) - (F} m (F> m )) ds, < I < K - 1, (2.16) 

21 



we see that 



K-l 



Wm,R < LC(R) + E Qr 



(2.17) 



so it remains to estimate E^^ ^ \ r)i\. We have 



\m\ 



< 
+ 
+ 



(F(F' m (s), <p u ' m (s)) - F{r< m {n), (p v ' m (ri) + z/^ 1 A m (s - n)))ds 
P 1 (Fil^in),^^ + I /- 1 A m (s - 77)) - (F} m (F' m ( n ))) ds 

I ((F) m (F' m (r ; ) - (F) m (F' m (s)) ds 

J T, 



--: T] + Tf + Tf . 



To estimate the quantities we first optimise the choice of the phase r . Consider 
the events £j, 1 < Z < K, 



Si = {F k {n) < 7 }, where 7 > a = 1/10. 



(2.18) 



By Lemma 2.1 and the Fubini theorem we can choose r G [0, L) in such a way 
that 



K- 1 ^P(^) = x( 7 - 1 ; J R,m). 



z=o 



For any I consider the event 



Q l = { sup \nr)-nn)\ hl > Pi(R)L 1/3 }, 

n<T<T l+1 

where Pi{R) is a suitable polynomial. It is not hard to verify (cf. [KP08]) that 
P(Qi) < XociL- 1 ). Setting 

Ti = Ei U Qi 



we have that 



K-l 



-1/2. 



m =: x. 



Accordingly, 



1 £ |(E - E Tl )Tj I < ^ P(7i) < P(R)x := xx, j = 1, 2, 3. 

«=0 1=0 
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Similar, since for uj g Q R the integrand in (2.16) is < Q(R), then 



(2-19) 



If to i tt R U T h then for r G [r,, r J+ i] we have that J£(t,) > 7 - Pi(i?)L 1/3 > §7, 
if z/ is small. This relation and (2.4), (2.5), (2.6) imply that 

PfWIv^s) - (V u ' m (ri) + ^A m (s - n)\ > v a for some s G fr,r m ]} 

< x 00 (z/- 1 ;i?,m). 

Accordingly, 

(2.20) 



It is clear that 



E E ^u^Tfj < P(R)L^ = P(Ry/°. (2.21) 

So it remains to estimate the expectation ofJ^Tf. For any oj QrUFi abbreviate 

F(4>) = F(i v > m (ti),<pr> m (ti) + 1>), i> G T m , 
where in the r.h.s. tp is identified with the vector (■?/>, 0, . . . ) G T°°. We can write 



Tff as 



T 2 — 
L l — 



rn+i 

/ F(u- l A m (s-n))ds- (F) 
J n 











= L 






LJo 



Since the function F(ip) is analytic and the vector A is non- resonant, then by 
Lemma 1.1 Tf < Lx(z/~ 1 L; m, R, 7, A). Therefore 



(2.22) 



Now (2.17), (2.19) and (2.20)- (2.22) imply that 

SI" < x 0O ( J R) + C(R)m- 1/3 + x(z/~ a ; i2, m) + x( 7 " 1 ; i2, m) 
+CV a + P(R)v l ' & + x(z/- 1 / 2 ; m, i2, 7, A). 

Choosing first R large, then m large and next 7 small and v small in such a way 
that (2.15) and (2.18) hold, we make the r.h.s. arbitrarily small. This proves the 
lemma. 
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2.4 Joint distribution of actions and angles. 

Denote \t v a = V(I v (s),(p v (s)) = (I x ip) o V(u u {s)), where u v {s), < s < T, is 
a solution of (1.1) and (I u ,(p v ) is a solution of the system (2.3), (2.4). For any 
/ G Li(0,T), / > 0, such that // = 1, set /2"(/) = £ f(s)jl u s ds. Also let us 
denote m°(/) = J* T f(s)V(I°(s)) ds; this is a measure on 

Theorem 2.4. For any / as above, 

ji u J(f) -± m°(f) x dy? as z/j ->■ 0. (2.23) 

Proof. For a piecewise constant function / the convergence follows from Theo- 
rem 2.2 and Lemma 2.3 since by the lemma, for any < T\ < T 2 < T, the integral 
jj 2 ^( /iy ( s )' ds is close to ff{F)(r(s)) ds, and by the theorem the inte- 

gral f£{F)(I v (s)) ds is close to jJ(F)(/°(s)) ds = f JOO F(I°(s),i/>) dsd^j (we 
are applying the lemma and the theorem on segments [0,T{\ and [0,T 2 ]). 

To get the convergence for a general function / we approximate it by piecewise 
constant functions. See Section 2 of [KuklO] for details. □. 

2.5 Multidimensional case. 

Let (2.1) be not eq. (1.1), but eq. (1.14), written in the f-variables. Now we should 
consider (2.1) as an equation in a space h r , r > d/2. The maps P l : h r — » h r and 
P 2 : h r ->■ /i r are smooth and the differentials gT'P 1 (v) : h r x ■ ■ ■ x h r h r are 
poly-linear mappings such that their norms are bounded by polynomials of \v |^r. 
This allows to apply to eq. (2.1) the methods of [KP08] 1 in the same way as in 
Sections 2.3-2.4 and establish validity of Theorems 2.2 and 2.4. 

3 Effective equations and uniqueness of limit. 

Let (2.1) be eq. (1.1) or eq. (1.14), written in the f-variables, and (2.11) - the 
corresponding averaged equation. Accordingly, by h we denote either the space 
h 1 as in Section 1, or the space h r , r > d/2, as in Section 1.3. For simplicity we 
assume that p and q in (1.14) are integers. If they are not, then in the calculations 
below the nonlinearities |w| 2p -a and |-a| 2,? -a should be modified by Lipschitz terms 
which cause no extra difficulties. 

was assumed in [KP08] that the relevant maps and vector-fields are analytic. This ana- 
lyticity was imposed only for simplicity. Sufficiently high smoothness and polynomial estimates 
on the corresponding high order differentials arc sufficient for all construction of [KP08]. 
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3.1 Effective equations. 

Let us write the averaged drift (vk • Pk) and the averaged diffusion (Ak r ) in the 
form (2.8) and (2.9), respectively. Using (2.2) we write the term R'(v) in (2.8) as 

3 . 3 

K(v) = E / e- i6k P™(*ev) d6 =: J2 k > 1. 

m=l m=l 

By (1.22) and (1.23), 

R\v) = x J (J^G^ev)^ d9 

= -xj <$>_ B A<$> B v d6 + x J <f>_e£ (<S> e v)d9 (3.1) 

= -xA + xR°(v), R°(v) = J ' $_gC°($gv)d6, 

since A commutes with the rotations The operator R° is bounded and self- 
adjoint in h°. For any v we have 

(R\v), v) = x J (AG&gv, G§ e v) d6 < -Cx\v\ 2 hl , (3.2) 

since ||G$^||i ~ \® e v\ h i = \v\ h i- Writing in (3.1) C°(v) = ^ oV o G(v) as 

C°(v) = xV(h 2 o G)(v), h\u) = ^J V(x)\u(x)\ 2 dx, 

we have R°(v) = V(h 2 o G)(v). Since 

(h 2 o G)(v) = \ J] / (V(x)e»>v m (x), e**vm{x)) dd 

= ]-^2\ Vl \ 2 M h M t = (Vw,^), 



i 

then R° = diag{M ; , I > 1}. Accordingly, 

R 1 = xdiag{-A^ + M u I > 1} > 0, M l = (Vtpi, <p t ) (3.3) 
The term R 2 is defined as an integral with the integrand 

$- e P 2 Mv) = -TflS-etffluM |«=Go* fl «=: ^(u). 
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Writing f p {\u\ 2 )u 2 as Vh p {u), where h p {u) = j F p {\u\ 2 )dx, (F p )' = and 
denoting G o $ e = L e , we have 



So 



F«,(u) = - lR L* e VhP{u) \ u=Le{v) = - lR V{h p oL e {v)). 

R\v) = - lR V v ( [ (h p o G)($ e v) ffl) = lR V v (K p o G). (3.4) 

Similar R 3 (v) = —i^iV v {h q oG) (since the operator Go$ e is complex-linear). 
As (h q o G) is a function solely of the actions (Ji, J 2 , • • ■ ), then V Vk (h q o G) G C 
is a vector, real-proportional to v k . Therefore v k ■ R\{v) = for each k. That is, 



((v k ■ P k ))(v) = v k ■ R\(y) + v k ■ R 2 k (v), 
where R 1 and R 2 are defined by (3.3) and (3.4). Now we set 

R(v) = R 1 (v) + R 2 (v) 
and consider the following system of stochastic equations: 

dv k (r) = R k (v)dr + Y k d(3 k , k > 1. 
Equations (3.6) are called the system of effective equations. 
Example 3.1. (p = 1). Now h\u) = \J \u\ A dx. So 



(3.5) 



(3.6) 



h 1 oG(v)= 1 -J 



^2v k ip k (x) 



Since 



(v kl v k2 v ks v k4 



IVki] \Vk 2 





dx=^ ^2 v kiVk 2 v k3 v k4 J (p kl (p k2 (p k3 tp k4 dx. 

ki,k2,k 3 ,k4 

\ 2 if k\ = k 3 , k<i = k 4 or k\ = k±, k<i = k 3 , 
otherwise, 



then 



ki^k-2 



2 if d = 1 and p is an integer, than / p (|w| 2 ) 



\2p 
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where L 'k lk2 = WIA dx - So that 

Rl(v) = ~lR^v k (h 1 o G)(v) = - lR v k ( \v k \ 2 L' kk + 2 \vi\ 2 L' k 

V l+k 

= -lRVk^2\vi\ 2 L M . 

I 

Here L kk = L' kk and L k i = 2L' kl if k ^ I. So the system of effective equations 
becomes 

dv k = -v k (x(X k -M k )+ lR J2\vi\ 2 L k i)dr + Y k d/3 k , k > 1. (3.8) 

i 

□ 

If v(r) = {ffc(r), A; > 1} satisfies (3.6), then for I k = I{v k {r)) we have 

dl k (r) =v k - R k (v)dr + Y k 2 dr + Y k v k -d(3 k , k > 1. (3.9) 

By (3.5) the drift in this system equals ((ffc • Pk)(I) + Y k )dr, while the diffusion 
matrix is ^rlffcl 2 ^ 2 = (A fcr ). So system (3.9) has the same set of weak (= 
martingale) solutions as (2.11), see [Yor74]. We have got 

Proposition 3.2. Let v{r) be a weak solution of (3.6) such that v(0) = vq and 
E sup \v(T)\ 2 h n < C\v \l n + C(n,T), Vn. (3.10) 

0<t<T 

Then Il/(t>(r)) is a weak solution of the system (2.11), satisfying (2.12) and such 
that 1(0) = I . 

That is, the solutions of eq. (2.11) which can be obtained as limits (when 
v — > 0) of actions I v (u(t)) of solutions for (1.1) (or (1.14)) are those which can 
be covered by "regular" solutions of (3.6). 

The 'right' inverse statement to Proposition 3.2 is given by the following 

Proposition 3.3. Let I°(t) be a weak solution of the averaged equations (2.11), 
constructed in Theorem 2.2. Then there exists a weak solution v°(t) of (3.6) such 
that v(0) = v , satisfying (3.10), and such that £>(n/(t; (-))) = V(I°{-)). 

For a proof we refer to Section 3 of [KuklO], where the assertion is established 
in a similar but more complicated situation. 
System (3.6) is invariant under rotations $e: 
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Proposition 3.4. Let v(r) be a weak solution of (3.6), satisfying (3.10). Then, 
for any 9 G T°° ; & e v(r) is a weak solution of (3.6), satisfying (3.10). 

Proof. Applying to (3.6) we get that 

d($ v) = $ e R(v) dr + <S> e Y d(3(r), Y = diag{F fe }. 

The vector fields i? 1 (f ) and R 2 (v ) both are obtained by averaging and have the 
form R j (v) = J &^ e F j (& e v) dO. So they commute with the rotations, as well as 
their sum R(v), and we have 

d($ e v) = R($ev) d6 + Yd($ e P(T)). 

Since P$g/3(r) = V(3{r), then the assertion follows. □ 



3.2 The uniqueness. 

Let v 1 (t) and v 2 (r) be solutions of the effective system (3.6). Denoting v = v 1 —v 2 , 
we have that 

l£\<T)\h < -Ml + W) - i?V) V - V 2 )- 

Consider the last term, denoting v 3 = Qgyi , u\ = G(v J e ). Since R 2 {y) is an 
integral over T°° with the integrand —^ R ^_ e ^(\ug\ 2p u e ), where u e = G(^g(v)), 
then 

{R 2 (v 1 )-R 2 (v 2 ),v 1 -v 2 ) = ->y R J (^(lul^ul-lu^ull^-^dB 

= -1R J ({\uTu) - \u 2 \^ul ul - u 2 ) dO. 

The integrand in the r.h.s. is non- negative. So 

1 d 



2dr 



v\io<-Cx\v\ 2 hl , (3.11) 



(i.e., the effective system (3.6) is strongly monotone). Therefore a strong solution 
of the system (3.6) is unique. By the Yamada-Watanabe argument (see [KS91]) 
a weak solution also is unique. We have got 
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Theorem 3.5. Let I u (t) = I(u u (t)), where u u (r) is a solution of eq. (1.1) or of 
eq. (1.14) and u u (0) = u . Then 

V{r{-)) ^ Q° as v 

m #ie space "H/, where Q° is a weak solution of (2.11), satisfying (2.12), (2.13). 
There exists a unique weak solution v{t) of the effective equations (3.6), satisfying 
(3.10), such that v(0) = V(uq) and £>(II 7 (>(-)) = Q° ■ 

4 Stationary solutions. 
4.1 Averaging. 

Again, let (2.1) be eq. (1.1) or eq. (1.14), written in the f-variables, and (2.11) 
be the corresponding averaged equation. Accordingly, by h we denote either the 
space h 1 as in Section 1, or the space h r , r > d/2, as in Section 1.3. Assume 
that the corresponding w-equation is regular in the space W (e.g., d — 1 or the 
assumptions, given at the end of the Section 1.3 are fulfilled), and that it has a 
unique stationary measure [i v (see Sections 1.2, 1.3). 

Let u' v [f) be a stationary in time solution of equation (1.1), T>(u' u (t)) = fi u . 
By estimates in Section 1 the set of laws V(I' U (-)), where l' v = I(u' u (t)), is 
compact in hj. Let Q' be any limiting measure as Vj — > 0. Clearly it is stationary 
in t. The same argument that was used to prove Theorem 2.2 (cf. [KP08]) imply 
that Q' is a stationary solution of the averaged equation: 

Proposition 4.1. The measure Q' is the law of a process I'ij), < r < T, which 
is a stationary weak solution of the averaged equation (2.11). It satisfies estimates 
(2.12), (2.13), and the stationary measure n = V(I'(Q)) meets estimates (1.11), 
(1.12) with v = 0. 

The measures (/ x tp) o fi u = V(I' I/ (s),ip" / (s)) satisfies (2.23) for the same 
reason as in Section 2.4. Since the measure [i v is independent from s, then now 

V(I' v (s),<p' v (s)) ^ 7T x dip as vj->0. (4.1) 

In the stationary case relation (2.7) implies that 

P{/'fc(r) < 5} ->• as 5^0, (4.2) 

uniformly in v. In particular, 

7r {/|4 = 0} = \fk. (4.3) 
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4.2 Lifting to effective equations. 

To study the limiting measure n further we lift it to a stationary measure of the 
effective system (3.6). We start with 

Lemma 4.2. Assume (1.10). Then the system (3.6) has a unique stationary 
measure m. 

Proof. Relation (1.10) implies that Yj, ^ for all k. That is, the noise in 
the effective equations is non-degenerate. Moreover, the coefficients Y k satisfy 
(2.10). Since solutions of (3.6) satisfy estimates (3.10) and since any two solutions 
converge by (3.8), then the assertion follows. E.g., see [KS10], Section 3. □ 

Let v(t) be a stationary solution of (3.6), T>(v(t)) = m. By Proposition 3.4, 
$o(v(t)) also is a (weak) stationary solution. So V($qv(t)) = §g o m is a sta- 
tionary measure for (3.6). Since it is unique, then 

$ e om = m V#GT°°. 

Accordingly, o m is a rotation-invariant measure on T°°, i.e. LT^ o m = dip. 
This implies that in the (I, (/?)-variables the measure m has the form 

dm = m 1 (dl) x dip (4.4) 

Proposition 4.1 applies for any time-interval [0, T]. So, replacing the sequence 
Vj — > by a suitable subsequence vy — > we construct a stationary process 
I'{ T )-i T > 0, such that I'"*' (t) converges to I'ij) in distribution on any finite 
time-interval. Using Theorem 3.4 and discussion at the end of Section 3.1 we 
construct a solution v'(t) of (3.6) such that V(TIj(v'(t)) = n. Since V(v'(t)) — 1 m 
as r — )> oo, then 

7r = ri/om. (4.5) 
That is, the measure n is independent from the sequence Vj. We have got 

Theorem 4.3. If (1.10) holds, then Iofi u ^n = H I om, where m is the unique 
stationary measure of the effective system. 

In view of (4.1), (4.4) and (4.5), 

(I x ip) o fi u — ^ (11/ x 11^) o m as v — > 0. 

Denote /i+ = G /i | v s ^ Vj}. By (4.2) and (4.3) (* o /i !/ )(/i+) = 1 and 
m(h + ) = 1. So the convergence above implies that 
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Theorem 4.4. // (1.10) holds, then fi u — ^ G o m as v — >■ 0. 

Example 4.5. (Hamiltonian perturbations.) If in (1.14) 7^ = 0, i.e. if the non- 
linear term of the perturbation is Hamiltonian, then the effective system is the 
linear equation 

dv(r) = R}(v) dr + Yd/3, 

where R l is defined in (3.3) and Y = diag{Y"fc, k > 1}. Let v(0) = 0. Then v(r) 
is the diagonal complex Gaussian process 

v(t) = [ T e (T - s ^ Rl Yd(3(s), R 1 = x(A - R°). 
Jo 

So the stationary measure for the effective system, Vv(oo), is a direct sum of 
independent complex Gaussian measures with zero mean and the dispersions 
x- 1 F fc 2 /(A fc -M fc ), k> 1. 

The fact that a Hamiltonian nonlinearity produces no effect in the first order 
averaging (i.e. for the slow time r < 1) is well known in the theory of weak 
turbulence. To produce a non-trivial effect, the Hamiltonian term — i'~fjf q (\u\ 2 )u 
should be scaled by the additional factor z/" 1//2 , and for the weak turbulence 
theory to apply to calculate this effect we should send the size of the x-torus to 
infinity when v — > 0, see [Nazll]. 

Example 4.6. (p = 1, continuation). If p = 1, then the effective equations become 

dv k = -v k (x(A fe - M k ) +1rJ2 M%z) dr + Y k d(3 k . (4.6) 

1 

Assume that the random force in (1.1) (or in (1.14)) is small and is mostly 
concentrated at a frequency j*. That is, 

bj,=e<l, 0<fy<£ if l^j*. 

Then the numbers Y k are of order e and are concentrated close to j*, i.e., 

Yj.~e, Y^sCnII-j^ VI, N. 

So if v(t) is a stationary solution of the effective equations and E k = |E|t>fc(r)| 2 , 
then 

E jm ~ e 2 \^\ E l < e 2 C N \^\l - j*\~ N V/, N. 

That is, the systems (1.1) and (1.14) exhibit no inverse or direct cascade of energy. 
For other polynomial systems (1.1) and (1.14) situation is the same. Certainly 
this is not surprising since by imposing the non-resonance condition we removed 
from the system resonances, responsible for the two energy cascades. 
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5 Equations with non-viscous damping. 



Following Debussche-Odasso [DO05] we now discuss equations (1.1) with non- 
viscous damping, i.e. with h = but with 7^ > and p = (Debussche-Odasso 
considered the case p — 0, q — 1): 

ii + iv' 1 (-u xx + V(x)u) = -j R u - i^u^u + -^^bjP^e^x), ^ 

u(x) = u(x + 2n) = —u(—x); 

u(0) = m . (5.2) 

Estimates (1.4), (1.8) and (1.9) are valid with x = 0. Jointly with an analogy of 
estimate (1.17) with x = 0, m = 1 they imply that for u G V, 2 the set of actions 
J"(r) = i"(^(r)) of solutions for (5.1) , (5.2) is tight in Uj. As in Section 2, 
any limiting measure Q° = \imT>(I l/j (-)) is a law of a weak solution I°(r) of the 
averaged equations (2.11)^ =0 with 1(0) — I — I(uq). Constructions of Section 3 
remain true, so I°(t) may be lifted to a weak solution v°(t) of the effective 
equations (3.6) x=0iP =o- Now R l = and, repeating constructions of Example 3.1 
we see that Rl(v) = —^rVu- So the effective equations become the linear system 

dv k (r) = -~f R v k dr + Y k d(3 k . (5.3) 

This system has a unique solution v(t) such that v(0) — vq — ^(uo). So 

]imV(F(.))=VU I (v(-)). 

Due to results of [DO05], eq. (5.1) has a unique stationary measure \i v '. 
Repeating arguments from Example 4.5, we see that when v — > 0, the measures 
^ o ji v converge to the unique stationary measure of eq. (5.3) which is 

r 



m = V f dmg{e~ slR Y k }d(3 k (s). 



This is a direct sum of independent complex gaussian measures with zero mean 
and the dispersion Y k 2 /^ R) k > 1. So every solution u(t) of (5.1) satisfies the 
Gaussian limit 

lim lim Du(r) = G o m. 

If we replace in (5.1) the linear damping by the nonlinear term — r yn\u\ 2 u, then 
the effective system (5.3) should be replaced by the nonlinear system (4.6) with 
\ k = M k = 0. In this case the limiting measure is non-Gaussian. 
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